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The finite field calculations are performed for two heavy frequency-standard candidates In + and 
Sr. The progressive hierarchy of electron correlations is implemented by the relativistic coupled- 
cluster and configuration interaction methods combined with basis set of increasing size. The dipole 
polarizabilities, dipole hyperpolarizabilities, quadrupole moments, and quadrupole polarizabilities 
are recommended for the ground state 5s 2 1 5'o and low-lying states 5s5p 3 .Po i 2 of In + and Sr. Com¬ 
parative study of the fully and scalar relativistic electron correlation calculations reveals the effect 
of the spin-orbit interaction on the dipole polarizabilities of In + and Sr. Finally, the blackbody- 
radiation shifts due to the dipole polarizability, dipole hyperpolarizability, and quadrupole polariz¬ 
ability are evaluated for the clock transition 5s 2 1 So - 5s5p 3 Pq of In + and Sr. 

PACS numbers: 31.15.ap, 31.15.aj, 32.10.Dk 


I. INTRODUCTION 

Polarizabilities describe the ability of atoms, ions, and 
molecules to be polarized in an external electric field, 
which is very useful in many areas of atomic and molec¬ 
ular physics M- In recent years, optical lattice clocks 
have shown a tremendous progress [5H8| , which has stim¬ 
ulated a great deal of interest in performing precision 
calculations of atomic polarizabilities and hyperpolar¬ 
izabilities. For instance, the difference of the static 
dipole polarizabilities between two associated states of 
an atomic clock transition determines the blackbody- 
radiation (BBR) shift that is crucial in evaluating the 
error of an atomic clock, which has been calculated for 
A1+ Cab Sr, Yb, In+, Sr + , Hg + , Mg, Ca and so on 
[3, y, 19M12II . The quadrupole polarizability is another im¬ 
portant quantity that is related to the multipolar BBR 
shift, which has been evaluated for Sr, Ca + , Sr + , and so 
on [ly, ll4Jj . As a higher-order response to the applied 
electric field, the hyperpolarizability also contributes to 
the energy shift of an optical frequency standard, being 
small but not necessarily negligible, which has been in¬ 
vestigated by Ovsiannikov, et al [jj|, [l|| . 

Among the current atomic optical clock candidates, an 
important category is based on the ns 2 1 Sq - nsnp 3 Pg 
optical transition [lf-yl, where the upper state is one of 
the metastablc triplet states 3 PJ with J=0, 1, and 2. 
Two such examples are the In + and Sr optical clocks us¬ 
ing the 5s 2 1 S'o - 5s5p 3 Pq transition @,[l(J. The polariz¬ 
abilities of such systems have been targeted with increas¬ 
ing theoretical efforts to perform precision calculations. 
For example, Sahoo et al. have calculated the ground- 
state dipole polarizability of Sr by using relativistic cou¬ 
pled cluster (CC) method H3- ' Safronova et al. have 
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given high accuracy data for the dipole polarizabilities of 
5s 2 1 S'o and 5s5p 3 Pq of In + and Sr using the configura¬ 
tion interaction (Cl) -fall-order method floL fl2l ] . Mitroy 
et al. have calculated the 5s 2 1 S’o and 5s5p 3 P° states 
of Sr by using Cl with a semiempirical core polarization 
potential [H, [l9| . Porsev et al. have calculated the Sr 
5s 2 1 So and 5s5p 3 Pq 1 states using the Cl method with 
many-body perturbation theory [20(. Very recently they 
have calculated the 5s5p 3 P° state by applying the more 
accurate Cl + all-order method [2l|. From these works, 
one can see that the most currently available dipole po¬ 
larizability data for In + and Sr are concentrated on the 
3 P 0 ° component, whereas the data for the 3 P° and 3 P^ 
components remain scarce. Besides, the quadrupole mo¬ 
ment of Sr in the 5s5p 3 P° state and its quadrupole po¬ 
larizabilities for the 5s 2 1 S and 5s5p 3 P° states have 
been calculated in earlier works by Mitroy, et al. us¬ 
ing the CICP method Emu and Porsev et al. using 
CI+MBPT and Cl+all-order methods [2l|, [13- However, 
the recommended data for the quadrupole moments and 
quadrupole polarizabilities of In + remain not available. 

Among various theoretical approaches in calculating 
polarizabilities, the finite-field method can provide reli¬ 
able data if the field-dependent energies are evaluated to 
a high precision. This method is often implemented in 
computational codes for atomic and molecular property 
calculations [23-l26j. In this method the employed ex¬ 
ternal field breaks the degeneracy of multiple states and 
thus the J and Mj-resolved polarizabilities of the nsnp 
3 Pj metastable states can be obtained directly without 
using basic vector algebra, where J and Mj are the total 
ang ular momentum and its magnetic quantum number 
[271 1 • I n particular, the influence of the spin-orbit cou¬ 
pling interaction on the J-resolved polarizability can be 
revealed through comparative studies between full and 
scalar relativistic calculations. Successful application of 
the finite-field method to atoms and ions of optical fre¬ 
quency standard has been demonstrated for Al + [28|, [29[ . 
However, the finite-field method has not been widely ap- 
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plied to calculate the polarizabilities of heavy atoms and 
ions, for which relativistic effects and electron correla¬ 
tions are expected to be significant. It remains challeng¬ 
ing in treating these effects in this finite-field approach. 

In this work, the dipole polarizabilities, dipole hy¬ 
perpolarizabilities, quadrupole moments, and quadrupole 
polarizabilities for the ground state 5s 2 1 S'o and the low- 
lying excited states 5s5p 3 Pg -l 2 of In + and Sr are calcu¬ 
lated by applying the finite-field method. The convergent 
hierarchies of electron correlation and the convergent ba¬ 
sis sets are adopted in the relativistic CC and Cl cal¬ 
culations in order to obtain properties of high accuracy. 
The dipole polarizabilities of the ground-state 5s 2 1 S 0 
and the excited state 5s5p 3 Pg , as well as the scalar and 
tensor polarizabilities of the excited states 5s5p 3 P° 2 , are 
obtained for In + and Sr. Comparisons with previously 
reported data for 5s 2 1 S'o and 5s5p 3 Pg show good agree¬ 
ment. The effect of the spin-orbit coupling interaction 
on the studied property is analyzed through comparing 
the J and L-resolved values of the 5s5p 3 Pq 12 states, 
where L is the total orbital angular momentum. Besides, 
we recommend the values of the dipole hyperpolarizabili- 
ties, the quadrupole moments, and the quadrupole polar¬ 
izabilities of In+ 5s 2 1 S'o and 5s5p 3 P§i 2 states and Sr 5s 2 
1 S'o and 5s5p 3 Pq states, for which the previous available 
data are rather scarce. Finally, the BBR shifts for the 
dipole polarizability, hyperpolarizability, and quadrupole 
polarizability are evaluated for the clock transition 5s 2 
1 Sq - 5s5p 3 P 0 ° of In+ and Sr. 


II. THEORY 

The energy shift of an atom or ion in a homogeneous 
electric held can be expressed as [3Qj 

A E d (F z ) = -aFU 2 - 7 P z 4 /24 - • • • , (1) 

where F z is the electric held strength along the z direc¬ 
tion, and a and 7 are, respectively, the dipole polarizabil¬ 
ity and dipole hyperpolarizability. The J-resolved scalar 
and tensor polarizabilities are given by 

^=2JTT£ 0(J ’ M ^ (2) 


Q J a =Q(J,\Mj\ = J)-Q j , (3) 

where Q denotes either a or 7 , and Q( J, Mj) denotes 
either a or 7 for each Mj component with Mj being the 
projection of the angular momentum J in the 0 direction. 
The L-resolved scalar and tensor polarizabilities are also 
dehned by Eqs. © and © with J and Mj replaced by 
L and Ml respectively. If the spin-orbit coupling inter¬ 
action is not strong so that it can be treated as a pertur¬ 
bation, the relationship between the J- and L-resolved 


polarizabilities can be established by using the LS cou¬ 
pling approximation [27] 

-J -L 

a = a , 

a a( 3 P°) = -a J a ?P2)/2 = -a L j2. (4) 

The actual relationship between the J- and L-resolved 
polarizabilities can deviate from Eq. © when the spin- 
orbit interaction is strong and the LS coupling scheme 
fails to describe the polarizabilities. 

In a pure quadrupole electric field, the corresponding 
energy shift is 

A E d {F zz ) = -~8F zz /2 - a 2 F 2 zz /8 - • • • , (5) 

where F zz is the electric field gradient in the z direction, 
and 8 and a 2 are the quadrupole moment and quadrupole 
polarizability, respectively. The quadrupole moment 9 
is calculated under the condition \Mj\ = J or \Ml\ = 
L. For the 5s5p 3 P state of In + and Sr, the L-resolved 
quadrupole moment is given for the state of 3 P° with 
Ml = 1 , and the J-resolved quadrupole moment is given 
for the state of 3 P 2 with Mj = 2. 

III. METHOD OF CALCULATION 

The electric field-dependent energy is calculated at 
different levels of theories. Dirac-Hartree-Fock (DHF) 
calculations are performed that generate the reference 
states as well as optimized atomic orbitals. These cal¬ 
culations are implemented by the SCF module in the 
DIRAC package [31[ . Both the Dirac-Coulomb Hamilto¬ 
nian and the Dyall Hamiltonian ||33j are employed. The 
Dirac-Coulomb Hamiltonian gives a complete description 
of the relativistic effects. We carry out the relativistic CC 
and Cl calculations that are built on the Dirac-Coulomb 
Hamiltonian and yield fully relativistic results. The rel¬ 
ativistic CC calculations are implemented by applying 
the MRCC suite [12] , and the relativistic Cl calculations 
are implemented by applying the the KRCI module in 
the DIRAC package. The Dyall Hamiltonian contains 
only the spin-free terms with the spin-orbit interaction 
neglected. We carry out the scalar relativistic Cl calcu¬ 
lations that are built on the Dyall Hamiltonian and yield 
the scalar relativistic results. The scalar relativistic Cl 
calculations are implemented by applying the LUCITA 
module in the DIRAC package. 

Due to escalating computational time for adopting 
higher level correlation methods, the generalized active 
spaces (GAS) technique [HJ is used to restrict the num¬ 
ber of correlated electrons and atomic orbitals, which 
makes the computation tractable. In the GAS, the Dirac- 
Fock orbitals are divided into the inner-core, outer-core, 
valence, and virtual shells. No excitation is allowed in 
the inner-core shells, but the excitation in the outer-core 
and valence shells can be defined to any order. In or¬ 
der to obtain convergent results, we choose the outer- 
core that includes increasing number of core electrons. 
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In particular, the 4 d, 4s4p4d, and 3d4s4p4d shells that 
consist of the ‘(corelO)’, ‘(corel8)’, and ‘core(28)’ outer- 
cores, respectively, are considered for In + . Similarly the 
4s4p, 3d4s4p, and 3s3p3d4s4p shells that consist of the 
‘(core8)’, ‘(corel8)’, and ‘(core26)’ outer-cores, respec¬ 
tively, are considered for Sr. For both In + and Sr, the 
valence shells are comprised of 5s5p. The virtual orbitals 
with the orbital energy larger than a given cut-off value 
are neglected in the correlation calculations. We have 
carried out the calculations with the different cut-off val¬ 
ues in order to check the convergence of the obtained 
results with the truncation of the virtual orbitals. 

In the following, a CC calculation is referred to as 
‘(core n)SD’ and ‘(core n)SDT’ that include single and 
double (SD) excitations and single, double, and triple 
(SDT) excitations from the outer-core and valence shells 
into the virtual orbitals, where n means the number of 
outer-core electrons that are involved in electronic cor¬ 
relation calculations. A Cl calculation is referred to 
as ‘(core n)SD(2in4)SDT’ and ‘(core n)SD(2in4)SDTQ’, 
where the outer-core shells are restricted to single and 
double (SD) excitations, and the valence shells, where 
two electrons are distributed in the four 5s5p orbitals, 
denoted by ‘(2in4)’, add excited electrons from the outer 
core, are restricted to SDT and single, double, triple, and 
quadruple (SDTQ) excitations into the virtual orbitals. 

The Dyall’s uncontracted correlation-consistent 
double-, triple-, and quadruple-^ basis sets are used, 
which are called XC, with X=2 , 3, and 4, respectively 
[HI, 36|. Each shell is augmented by two additional 
diffused functions. The exponential coefficients of the 
augmented functions are determined according to 



where (n and (qv-i are the two most diffused exponents 
for the respective atomic shells in the original basis set. 
Arbitrary 4-8 finite field strengths are chosen in the range 
of F z =(0, 4.5xl0 -3 ) and of F zz =( 0, 4.5xl0 -5 ) in atomic 
units. The fitting is checked to remove the dependence of 
the properties studied on sampling. In our calculations, 
the criterion for energy convergence is set to be 10 -1 ° 
Hartree. 

We use the composite scheme (2g| to give the final value 
TFinai of a studied property. In the CC calculation, 

Pf inal = PSD + AP T + AP core , (7) 

where Psd is the value calculated for property P using 
the (corelO)SD method with the X=A(, basis set, A Pt 
and AP C ore are the corrections due to the triple excitation 
and more outer-core electrons. In the Cl calculation, 

PFinal = PsDT + A Pq + A P core , (8) 

where Psdt is the value calculated for property P using 
the (corelO)SD(2in4)SDT method with the X=A( basis 
set, and A Pq is the correction due to the quadruple exci¬ 
tation. The virtual orbitals with the orbital energy larger 


than 20 a.u. for In + and 10 a.u. for Sr are neglected in 
the correlation calculations. We also carry out the cor¬ 
responding calculations of the cut-off value of the virtual 
orbitals of 100 a.u. for In + and Sr. The truncation of vir¬ 
tual orbitals of the cut-off value 100 a.u. results in change 
of a J less than 0.5% in comparison with the cases of the 
cut-off value 20 a.u. for In + and 10 a.u. for Sr. The con¬ 
tribution of the virtual orbitals with the orbital energy 
larger than 20 a.u. for In + and 10 a.u. for Sr is therefore 
omitted in the composite scheme. 

The uncertainty in PFinal is mainly caused by three 
possible error sources. The first error is due to the finite 
basis set used for calculating Psd and Psdt ■ The con¬ 
vergence of P with respect to the basis set is very quick 
when the basis set is larger than X=3(), as shown by the 
A1+ results (29|. Thus, we assume empirically that the 
error for Psd and for Psdt is equal to half of the dif¬ 
ference of the P values calculated by using X=4£ and 
3£ basis sets. This error considers the possible correc¬ 
tion with respect to the value computed in the infinite 
basis. The second error comes from the estimation of 
A Pt and APq . Previous experience [29| has shown that 
the A Pt and APq contributions computed even with the 
smaller 2Q basis set are never in error by more than 50% 
with respect to the basis-set limit, and hence the error 
bar of APr and of APq is taken to be half of itself, i.e., 
0.5 x A Pt and 0.5 x APq . The third error is due to the es¬ 
timation of AP core . As we will show in the following, the 
studied properties start to converge even with medium- 
size outer cores such as core(18), which indicates that the 
error of AP core is not larger than P( core2 s) - -P( CO rei8) for 
In + and P( core 26 ) — -P(corei8) f° r Sr. Hence, the overall 
uncertainty in PFinal can be estimated to be Root-Mean- 
Square(RMS) of such three errors. 

Throughout this paper, atomic units (a.u.) are used, 
unless otherwise stated. The atomic units of a, a 2 , 7, 
F z . and F zz are, respectively, 1.648778 x 10 -41 C 2 m 2 J -1 , 
4.617048 x 10 _62 C 2 m 4 J _1 , 6.235378 x 10" 65 C 4 m 4 J- 3 , 
5.142250xl0 9 V/cm, and 5.142250xl0 n V/cm 2 . 

IV. RESULTS AND DISCUSSION 

A. Dipole polarizability and hyperpolarizability 

Tabled summarizes the results of a J for the 5s 2 1 Sq 
and 5s5p 3 Pq 12 states of In+ calculated by using the rel¬ 
ativistic CC method. Firstly, the obtained values of a J 
in the (corelO)SD calculation with the X=2(, 3£, and 4£ 
basis sets show a good convergence. The error in Psd 
is only about 0.01^0.06, which implies that the correc¬ 
tion with respect to the infinite basis set is very small. 
Then, the AP core correction is estimated with the X=3£ 
basis set. Upon inclusion of the 4s4p4d core electrons 
into the (corel8)SD calculation, a J decreases in compar¬ 
ison with the (corelO)SD case. The effect of adding more 
core electrons, 3d4s4p4d, is illustrated by the rather small 
difference between the (core28)SD and (core!8)SD ealeu- 
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TABLE I: Dipole polarizabilities a J of In + for the states of 5s 2 1 So and 5s5p 3 -Pop ,2 calculated by the relativistic CC method, 
where (corelO), (corel8), and (core28) correspond respectively to the 4 d, 4s4p4d, and 3d4s4p4d core shells included in the 
electron correlation calculations. 


Level of excitation a 

1 So 

3 po 

M) 

3 po 

\Mj\—0 

\Mj\=l 

3 po 

|Afj|=2 

a J 

Oil 

Basis:2C(23s,17p,13d,4f) 









(corelO)SD 

24.83 

27.91 

29.12 

32.34 

31.47 

28.88 

30.61 

-1.73 

(corelO)SDT 

24.54 

27.50 

28.75 

32.26 

31.13 

28.42 

30.27 

-1.85 

A P T 

-0.29 

-0.41 

-0.37 

-0.08 

-0.35 

-0.46 

-0.34 

-0.12 

Error in APt 

±0.15 

±0.21 

±0.19 

±0.04 

±0.18 

±0.23 

±0.17 

±0.06 

Basis:3C (30s, 23p, 17d, 5f ,3g) 









(corelO)SD 

24.90 

27.01 

27.98 

30.75 

30.09 

28.10 

29.42 

-1.33 

(corel8)SD 

24.71 

26.81 

27.83 

30.68 

30.09 

27.84 

29.31 

-1.47 

(core28)SD 

24.67 

26.77 

27.79 

30.63 

29.93 

29.79 

29.21 

-1.42 

APcore 

-0.24 

-0.25 

-0.19 

-0.12 

-0.16 

-0.30 

-0.21 

-0.09 

Error in A P CO re 

±0.04 

±0.04 

±0.04 

±0.05 

±0.16 

±0.05 

±0.09 

±0.05 

Basis:4C(35s,29p,20d,7f,5g,3h) 









(corelO)SD,P s.d 

24.86 

26.91 

27.87 

30.64 

29.98 

28.01 

29.32 

-1.31 

Error in Psd 

±0.03 

±0.05 

±0.06 

±0.06 

±0.05 

±0.04 

±0.05 

±0.01 



PFinai- 

= ■Ps D +A Pcore + A Pt 




Final data, PFinai 

24.33 

26.25 

27.31 

30.44 

29.48 

27.25 

28.78 

-1.53 

Uncertainty (%) 

0.62 

0.82 

0.73 

0.28 

0.82 

0.87 

0.69 

4.99 

Ref. [Toll 

24.01 

26.02 








a The relativistic Cl calculation is performed using the 3£ basis set at (corelO)SD(2in4)SD<2 level, which yields 
a:=25.06 and 27.93 for 1 So and 3 Pq , a J =28.59 and 0^=0.34 for 3 P°, and a J =30.30 and a^=—1.38 for 3 P^• 


TABLE II: Dipole hyperpolarizabilities r y J of In + for the states of 5s 2 1 S'o and 5s5p 3 P§,i ,2 calculated by the relativistic CC 
method, where (corelO), (corel8), and (core28) correspond respectively to the 4d, 4s4p4d, and 3d4s4p4d core shells included 
in the electron correlation calculations. 


Level of excitation a 

1 So 

3 po 

3po 

Mj | = 1 



3 po 

P2 



\Mj\=0 


|Mj|=2 

7 J 

7 a 

Basis:2C(23s,17p,13d,4f) 









(corelO)SD 

3695 

14752 

21119 

27116 

20729 

7294 

16632 

-9338 

(corelO)SDT 

3640 

15134 

21611 

27546 

20989 

7773 

17014 

-9241 

AP t 

-55 

381 

492 

431 

260 

479 

382 

97 

Error in APt 

±28 

±191 

±246 

±216 

±130 

±240 

±191 

±49 

Basis:3C(30s,23p,17d,5f,3g) 









(corelO)SD 

3143 

13435 

19265 

26644 

20951 

6820 

16437 

-9617 

(corel8)SD 

3072 

13467 

19368 

26472 

20755 

6968 

16384 

-9645 

(core28)SD 

3065 

13464 

19327 

26480 

20703 

7056 

16400 

-9344 

AP C ore 

-78 

28 

62 

-164 

-248 

236 

-38 

273 

Error in A P COTe 

±7 

±3 

±42 

±9 

±52 

±88 

±16 

±72 

Basis:4£(35s,29p,20d,7f,5g,3h) 









(corelO)SD, Psd 

3122 

13057 

19125 

26265 

20571 

6765 

16187 

-9422 

Error in Psd 

±n 

±189 

±70 

±190 

±190 

±28 

±125 

±97 



-^Final - 

= Psd APcore +APt 




Final data, PFinai 

2989 

13467 

19679 

26532 

20583 

7479 

16531 

-9052 

Uncertainty (%) 

1.01 

1.99 

1.32 

1.08 

1.15 

3.43 

1.38 

1.44 


a The relativistic Cl calculation is performed using the 3£ basis set at (corelO)SD(2in4)SD<2 level, which yields 
a:=2715 and 16164 for 1 S'o and 3 Pq , a J =17247 and cx. J a =4263 for 3 P°, and a ,/ = 18614 and o:^=—8518 for 3 P £. 


lations, indicating that a J has entered into convergence 
region. Thus, AP core is estimated as the difference be¬ 
tween the (core28)SD and (corelO)SD calculations, which 
is about 0.09^0.30. By comparison of the (corelO)SD, 
(core!8)SD, and (core28)SD calculations, we see that the 


minimum number of core electrons that needs to be cor¬ 
related for In + is 18; in other words, the 4s4p4d core 
electrons need to be correlated in order to obtain accu¬ 
rate a J . Next, APp is estimated by the difference be¬ 
tween the (corelO)SD and (corelO)SDT values with the 
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TABLE III: Dipole polarizabilities a L and hyperpolarizabilities , y L for the 5s 2 1 S and 5s5p 3 P° states of In + obtained by the 
scalar relativistic Cl method, where (corelO), (corel8), and (core28) correspond respectively to the 4d, 4s4p4d, and 3d4s4p4d 
core shells included in the electron correlation calculations. 


Level of excitation 



O' 





7 



1 S 0 


3 

po 





3 po 


| M l |=0 

|Mi|=l 

a L 


|M l |=0 

|Mi | — 

=1 7 L 

7 a 

Basis:2C(23s,17p,13d,4f) 











(corelO) SD(2in4)SDT 

24.42 

32.24 

27.98 

29.40 

-1.42 

2832 

35528 

7099 

16576 

-9476 

(corel8)SD(2in4)SDT 

24.28 

32.23 

27.84 

29.31 

-1.46 

2521 

35712 

6972 

16552 

-9580 

(core28)SD(2in4)SDT 

24.26 

32.19 

27.82 

29.28 

-1.46 

2567 

35679 

6967 

16538 

-9571 

(corelO) SDT(2in4)SDTQ 

24.41 

32.17 

27.82 

29.27 

-1.45 

3650 

35577 

7381 

16780 

-9398 

A Pq 

-0.015 

-0.07 

-0.16 

-0.13 

-0.03 

789 

49 

282 

204 

78 

Error in Pq 

±0.008 

±0.035 

±0.08 

±0.065 

0.015 

±395 

±25 

±141 

±102 

±39 

AP C ore 

-0.16 

-0.05 

-0.16 

-0.12 

-0.04 

-294 

151 

-132 

-38 

-94 

Error in P core 

±0.02 

±0.05 

±0.0.02 

±0.03 

±0.01 

±46 

±33 

±5 

±14 

±9 

Basis:3C(30s,23p,17d,5f,3g) 











(corelO)SD(2in4)SDT 

24.38 

30.43 

27.17 

28.26 

-1.08 

2182 

34283 

6742 

15922 

-9180 

Basis:4C(35s,29p,20d,7f,5g,3h) 











(corelO)SD(2in4)SDT, Psdt 

24.34 

30.35 

27.12 

28.20 

-1.08 

2270 

34257 

6761 

15926 

-9165 

Error in Psdt 

±0.02 

±0.04 

±0.03 

±0.03 

±0.005 

±44 

±13 

±9 

±2 

±8 





Pf i 

inal = Ps DT + A P CO re 

±APq 




Final data, PFinai 

24.16 

30.22 

26.80 

27.94 

-1.14 

2765 

34457 

6910 

16092 

-9182 

Uncertainty (%) 

0.13 

0.23 

0.33 

0.28 

1.56 

14.46 

0.13 

2.05 

0.64 

0.44 


X = 2£ basis set, which is around —0.08 ~ —0.46. Both 
the A Pt and AP cole corrections are large in comparison 
with the error in Psd , and thus an omission of such a 
correction would lead to an underestimation of a J . The 
reference data for the 5s 2 1 S'o and 5s5p 3 Pq states of In + 
are also given in Tabic |U as calculated by Safronova et al. 
by using the Cl+all-order method M- Our results are 
consistent with their data with a discrepancy around 1 %. 
Finally, we recommend that d‘ / =27.31 for 5s5p 3 P° with 
\Mj\=l and d J =28.78 and a J a = —1.53 for 5s5p 3 P 2 . 

Our calculated results for r ) J are listed in Table HD The 
7 J values show a similar convergence trend with respect 
to the size of basis set. The values of A Pt, and AP cole 
are also larger than the error in Psd for r ) J , similar to 
the case of a J . As shown in Tables Hand El the largest 
source of uncertainty in the final data comes from the 
error in A Pt for both a J and 7 U As mentioned above, 
A Pt has important contribution to the final data. In this 
work, due to very high computational demand in using 
larger basis sets, A P T is calculated only with X = 2£. 
This basis set is much smaller, which may cause an over¬ 
estimation of the uncertainty of A Pt, like the cases of a J a 
for 5s5p 3 P 2 and r y J for 5s5p 3 P 2 with Mj = 2. Using 
the same composite scheme of convergence, we arrive at 
the recommended values of 7 J =2989, 13467, and 19679 
for the 5s 2 1 5'o, 5s5p 3 Pq , and 5s5p 3 P°, \Mj\=l states, 
respectively, and y J =16531 and 7 ^ = —9052 for the 5s5p 
3 P 2 ) state. The uncertainties of the final values of 7 are 
slightly larger than those for a J , which is comprehensible 
because the hyperpolarizability is a higher-order response 
that is more sensitive to a small energy variation and thus 
various contributions bring substantial corrections to y -7 . 
The relativistic Cl calculation is implemented at the level 


of (corelO)SD(2in4)SD<2, where the cutoff for the vir¬ 
tual orbitals is 2 a.u. in energy, and the basis set is 3£, 
as indicated at the footnotes of Tables H and El 

The results of a L and y L for the 5s 2 1 S I and 5s5p 
3 P° states of In + , calculated by the scalar relativis¬ 
tic Cl method, are contained in Table mil Ps DT is 
determined using the X=4( basis set, and A Pq and 
AP core are determined using the X=2( basis set as 
the differences between the (corelO)SD(2in4)SDT 
and (corelO)SDT(2in4)SDTQ, and between the 
(corelO)SD(2in4)SDT and (core28)SD(2in4)SDT, 
respectively. The triple excitation is considered in 
the scalar relativistic Cl calculations, which reduces 
the uncertainties in the final values of a L and 7 ^, 
as compared with the relativistic CC calculations. In 
Table Ell the A Pq and AP core corrections are also larger 
than the error in Psdt- This trend is in accordance 
with that in the relativistic CC calculations. Because 
A Pq is computed only with a much smaller basis set, 
an estimation of error A Pq can not be exact (like A Pt 
in Tables Hand El), which leads to an anomalously large 
uncertainty in the final data for the 5s 2 1 S'o state. Up 
to now, we can find that, for In + , increasing basis set 
up to X = 4£ has decreased the error in Psd with 
respect to the infinite basis set. However, the A Pt and 
APcore corrections are larger than the error in Psd- In 
this situation, both APp and AP core become crucial for 
an accurate evaluation of the dipole polarizability and 
hyperpolarizability of In + . 

Table IIVI summarizes the results of a J for the states 
of 5s 2 1 S'o and 5s5p 3 Pq x 2 of Sr, calculated by using the 
relativistic CC method. The relativistic Cl calculations, 
implemented at the (core8)SD(2in4)SD<2 level with the 
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TABLE IV: Dipole polarizabilities a J for the states of 5s 2 1 S o and 5s5p 3 P§,i ,2 of Sr calculated by the relativistic CC method, 
where (core8), (corel8), and (core26) correspond respectively to the 4s4p, 3d4s4p, and 3s3p3d4s4p core shells included in the 
electron correlation calculations. 


Level of excitation a 

1 

3 po 

3po 



3 po 

P2 




M) 

|Mj|=l 

|Mj|=0 

|Mj| = l 

|Mj|=2 

a J 

Oil 

Basis: 2£(23s,17p,12d,3f) 
(core8)SD 

205.1 

349.3 

377.1 

418.7 

397.1 

331.7 

375.3 

-43.5 

(core8)SDT 

200.1 

348.8 

379.5 

422.2 

499.8 

332.0 

377.2 

-45.2 

A P T 

-5.0 

-0.5 

2.4 

3.5 

2.8 

0.2 

1.9 

- 1.7 

Error in APt 

± 2.5 

±0.3 

± 1.2 

± 1.8 

± 1.4 

± 0.1 

± 1.0 

± 0.9 

Basis:3C(31s,22p,15d,7f,3g) 

(core8)SD 

204.8 

423.8 

460.8 

515.1 

490.5 

415.7 

465.5 

- 49.8 

(corel8)SD 

204.6 

415.3 

450.3 

503.8 

478.8 

406.3 

454.8 

-48.5 

(core26)SD 

204.7 

418.8 

454.4 

509.0 

486.9 

410.1 

460.6 

-50.5 

APcore 

-0.2 

-5.0 

- 6.4 

- 6.2 

-7.6 

-5.6 

-6.5 

0.9 

Error in AP core 

± 0.01 

± 3.5 

± 4.0 

± 5.2 

±4.1 

± 3.9 

± 4.2 

± 0.3 

Basis:4C(35s,27p,17d,9f,7g,3h) 
(core8)SD, Psd 

204.8 

446.4 

484.9 

546.6 

520.0 

441.9 

494.1 

-52.2 

Error in Psd 

0.0 

± 11.3 

± 12.1 

±15.7 

±14.8 

±13.1 

±14.3 

± 1.2 



PFinal : 

= PSD+ APcore 

+APt 





Final data,Ppinai 

199.7 

444.1 

480.9 

543.9 

519.2 

436.5 

491.1 

- 54.5 

Uncertainty (%) 

1.2 

2.7 

2.7 

3.1 

3.0 

3.1 

3.1 

2.8 

rcc m 

199.7 








CI+MBPT [13 

197.2 

457.0 

498.8 






CI+all-order [T^| 

194.4 

441.9 







cicp m 

204.5 


497.0 (27.7) b 






CI+all-order [Hj 



459.2 (26.0) 6 






Expt. [331 



(24.5) b 







a The relativistic Cl calculation is performed using the 3£ basis set at (core8)SD(2in4)SD<2 level, which yields 
a=179.5 and 394.5 for 5s 2 1 So and 5s5p 3 Pq , a^=409.7 and a^—23 for 5s5p 3 P°, and d^=455.2 and a^=—51.8 
for 5s5p 3 P% • 

b Here given are the scalar and tensor (in parentheses) polarizabilities for 5s5p 3 P°. 


TABLE V: Dipole hyperpolarizabilities 7 J for the states of 
5s 2 1 S'o and 5s5p 3 Pq of Sr obtained by the relativistic CC 
method, where (core8), (corel8), and (core26) correspond re¬ 
spectively to the 4s4p, 3d4s4p, and 3s3p3d4s4p core shells 
included in the electron correlation calculations. 


Level of excitation 

1 So 

3 po 

M) 

Basis:2C(23s,17p,12d,3f) 



(core8)SD 

484500 

7050552 

(core8)SDT 

510144 

6905904 

AP t 

25644 

-144648 

Error in APt 

±12822 

±72324 

Basis:3C(31s,22p,15d,7f,3g) 



(core8)SD 

680965 

3568344 

(corel8)SD 

672720 

3452688 

(core26)SD 

674592 

3289040 

AP C ore 

-6373 

-279304 

Error in AP core 

±1872 

±163648 

Basis:4C(35s,27p,17d,9f,7g,3h) 



(core8)SD, Psd 

672686 

3652171 

Error in Psd 

±8279 

±83827 

PFinal = PsD+AP cor e+APT 


Final data, Ppinai 

691957 

3228219 

U ncert ainty (%) 

2.22 

6.12 


TABLE VI: Dipole polarizabilities a L for the states of 5s 2 
1 S and 5s5p 3 P° of Sr obtained by the scalar relativistic Cl 
method, where (core8), (corel8), and (core26) correspond re¬ 
spectively to the 4s4p, 3d4s4p, and 3s3p3d4s4p core shells 
included in the electron correlation calculations. 


Level of excitation 

1 S 0 


3 

po 


\m l \= 

=0 \Ml | = 

1 ol L 


Basis:2C(23s,17p,12d,3f) 






(core8)SD(2in4)SDT 

199.5 

443.7 

312.6 

356.3 

- 43.7 

(corel8)SD(2in4)SDT 

198.2 

427.5 

302.3 

344.0 

-41.8 

(core26)SD(2in4)SDT 

197.8 

427.6 

302.3 

344.1 

-41.7 

(core8)SDT(2in4)SDTQ 

198.7 

457.2 

325.2 

369.3 

-44.0 

AP q 

-0.7 

13.5 

12.6 

12.9 

-0.3 

Error in APq 

±0.4 

±6.8 

± 6.3 

± 6.5 

± 0.2 

APcore 

-1.6 

-16.2 

-10.3 

-12.3 

2.0 

Error in A P COTe 

± 0.35 

± 0.01 

±0.05 

±0.04 

±0.04 

Basis:3C(31s,22p,15d,7f,3g) 






(core8)SD(2in4)SDT 

198.4 

528.4 

381.9 

430.7 

-48.8 

Basis:4<^(35s,27p,17d,9f,7g,3h) 





(core8)SD(2in4)SDT, P S d 

197.1 

547.7 

396.7 

447.0 

-50.3 

Error in Psd 

±0.7 

±9.6 

±7.4 

± 8.1 

±0.7 



^Final =■ Ps D T + A P CO re + A Pq 

Final data, PFinai 

194.7 

544.9 

398.9 

447.6 

-48.7 

U ncert ainty (%) 

1.3 

1.8 

2.5 

2.1 

3.5 
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TABLE VII: Quadrupole moments 9 and polarizabilities 02 , 
where the L-resolved values are obtained by the scalar rela¬ 
tivistic Cl method and the J-resolved values are obtained by 
the relativistic CC method for 5s5p 3 Pg 2 and the relativistic 
Cl method for 5s5p 3 P°. 




e L 

e J 

In + 

1 S 




3 P° 

4.36 

4.64 ( 3 P°) 

Sr 





3 po 

15.56 

15.76 ( 3 P 2 °) 

15.6 [H, 15.52 m 



“2 


In+ 

x s 

127 

129 


3 po 

145 

1425 ( 3 P 0 °) 

1678 ( 3 P 1 °) 

-859.3 ( 3 P 2 °) 

Sr 

x s 

4688 

4608 



4640 El 

4545 [H 


3po 

6756 

9.75Xl0 4 ( 3 P 0 °) 



6949 |ll 

1.17X10 5 ( 3 P 1 °) 
1.05x10 s [H 
-7.39X10 4 ( 3 P“) 


3£ basis set, are presented at the footnote area for com¬ 
parison. As seen from the (core8)SD data with the X=2 , 
3, and 4£ basis sets in Table HV1 the changes of a J for 
Sr 5s 2 1 S'o state are very small , while the a J values for 
Sr 5s5p 3 P 0 °i 2 states tend to increase obviously for the 
larger basis sets, implying the basis set effect is rather 
large for these states. The error in Psd for the 5s5p 
3 Pq 12 states of Sr is about 12.1^15.7, except for a J a of 
5s5 p 3 P£, which is larger than A Pt and AP core . This 
indicates that the primary dominant factor that affects 
the final a J of Sr comes from the size of the basis set. 
This trend differs from the results mentioned earlier for 
In + . Note that the error in A Psd is very small and is 
far less than A Pt and AP core for the case of In + . The 
possible reason for such a difference is because of the dif¬ 
ferent electron distribution of In + and Sr. As a positively 
charged ion, in general In + has compact electron distri¬ 
bution around nucleus. However, as a neutral atom, Sr 
would have more diffused electron distribution than In + 
in an external electric field, indicating that larger basis 
sets are required for computing Sr. For Sr, the uncer¬ 
tainty in a J is around 1.2~3.1, most of which is from the 
error in Psd- This reflects that the accurate calculation 
of a of Sr replies greatly on the basis set. The values of 
7 J for the 5s 2 1 S'o and 5s5p 3 Pq states of Sr, calculated 
by using the relativistic CC method, are summarized in 
Table [V] As a higher-order property, the final data of r ) J 
for the 5s 2 1 S'o and 5s5p 3 Pg states of Sr show a large 
dependence on A Pt and AP core , and also show larger 
uncertainties than that for a J . 

The previously reported dipole polarizabilities of Sr are 
also listed in Table HVl for comparison. They include the 


calculation of Sahoo et al. using the relativistic couple 
cluster (RCC) method na , the calculation of Porsev et 
al. using the Cl method with many-body perturbation 
theory (CI+MBPT) [20|, the calculation of Safronova et 
al. using the Cl + all-order method [lU, the calculation 
of Mitroy et al. using the Cl method with semiempiri- 
cal core polarization potential (CICP) 0 , and the cal¬ 
culation of Porsev et al. using Cl + all-order method 
[inf . Noteworthy to mention are the CI+all-order re¬ 
sults of Safronova et al. that give a=194.4 and 441.9 
for 5s 2 1 Po and 5s5p 3 Pq , respectively. For such two 
states, we obtain a J =199.7 and 444.1 in the relativistic 
CC calculation, which are slightly larger than the values 
of Safronova et al. The possible reason for such discrep¬ 
ancies may be due to underestimation of the magnitudes 
of A Pt and A P core for which our basis sets are much 
smaller. There are fewer reported data available for the 
5s5p 3 P° state. Porsev et al. have obtained a J =498.8 
for the 5s5p 3 P° state with Mj = 1 state by using the 
CI+MBPT method, which is close to our calculated re¬ 
sult a J =480.9. The scalar and tensor polarizabilities of 
5s5p 3 P° have also been given as d J =497.0 and a^=27.7 
by Mitroy et al. [3 and d J =459.2 and a„=26.0 by Por¬ 
sev et al. (2lj| . The experimental value of a J a for 5s5p 
3 P° is 24.5. By using the relativistic Cl method we ob¬ 
tain d J =409.7 and a^=23.0 for Sr 5s5p 3 P°, which dif¬ 
fers from Porsev’s values by 11%. Such a relativistic Cl 
calculation is implemented with much smaller basis set 
and has about 10% error in a J in comparison with our 
relativistic CC calculation. Considering this error, our 
results are much closer to the values of Porsev et al. [U . 
There is no reported data available for the dipole polariz¬ 
ability of the Sr 5s5p 3 P^ state. In the present relativistic 
CC calculations we obtain d J =491.1 and a^=-54.5 for 
5s5p 3 P£- For the hyperpolarizabilities of Sr we recom¬ 
mend 7 J =691957 and 3228219 for 5s 2 1 S'o and 5s5p 3 Pq, 
respectively. 

The results of a L for the states of 5s 2 1 S I and 5s5p 
3 P° of Sr are given in Table fVTl In the scalar relativistic 
Cl calculations, a L for the 5s 2 1 S' ground state seems to 
be insensitive to the basis set, APq, and AP core ; on the 
other hand, a L for the 5s5p 3 P° state is sensitive to these 
factors. a L of the 5s5p 3 P° state increases with the size 
of basis set and it reaches convergence at X=4 £. A Pp 
and AP core are substantial in magnitude but they are in 
opposite sign, resulting in a cancellation with each other 
in the final value of a L . Notice that this phenomenon is 
also found in Table HVl for the relativistic CC calculation 
but the cancelation is not as strong as that in Table IVII 
for the scalar relativistic Cl calculation. Combining the 
relativistic CC and scalar relativistic Cl results, we can 
conclude that the basis set has very important effect on 
the final values of a J and a L . A further expansion of 
the basis set will no double improve the accuracy of a. 
As shown in Tables m and IVII for both the relativistic 
CC and the scalar relativistic Cl calculations, AP core has 
reached convergence upon inclusion of 3d4s4p into the 
electron correlations. This indicates that the minimum 








number of the core electrons is 18 for the case of Sr, 
i.e., the correlation contribution from 3d4s4p needs to be 
included in all-electron calculation in order to achieve an 
accurate dipole polarizability. 

The relativistic effect in the four-component relativis¬ 
tic formalism can be understood as a combination of spin- 
orbit coupling effect and contraction or decontraction of 
the radial electron density, i.e., the scalar relativistic ef¬ 
fect. Among the 5s5p 3 P° states, the 3 Pq component 
has a spherically symmetric electronic density. In this 
case, the effect of the pure spin-orbit interaction on the 
polarizability can be reflected by comparing the results 
obtained by fully and scalar relativistic calculations for 
the 3 Pq component, as we have done for Al + , in terms of 

Q L ( 3 P°) - Q J ( 3 PS ) 

WWW) ' (9) 

Such a fractional difference is about 8 % in a and 22% 
in 7 for the 5s5p 3 P° states of In + , showing an evident 
spin-orbit effect in a and 7 . The fractional difference 
[a L — a J ( 3 Pg)\/a J ( 3 Pq) amounts to 2% for Sr, which is 
smaller than that of In + . This implies that the effect of 
the spin-orbit coupling on a of Sr is weaker than In + . 
The relativistic effect can be elucidated by the fractional 
difference as follows, 


out with the AT=3C basis set. From Table IVUl we can 
found that, the values of 9 L and 9 J for the 5s5p 3 P° state, 
as well as the values of a 2 and a 2 for the 5s 2 1 S state, 
are in good accord, which indicates that these properties 
are insensitive to the spin-orbit coupling. In contrast, the 
values of a 2 and a 2 for the 5s5p 3 P° state show clear 
discrepancy, and a 2 changes greatly for different J com¬ 
ponents. This means that the spin-orbit coupling plays 
an important role on the quadrupole polarizabilities of 
the 5s5p 3 P° states of In+ and Sr. 

The quadrupole moments of Sr in the 5s5p 3 P° state 
were calculated in earlier works by using the CI+MBPT 
method |22j and the CICP method [23| . Our results agree 
with their data within 1%. The quadrupole polarizabili¬ 
ties for the 5s 2 X S and 5s5p 3 P° states of Sr were calcu¬ 
lated using the CICP method [l 8 j and the CI+all-order 
method with the random phase approximation (RPA) 
applied [21]. For the ground state 5s 2 1 S, our results are 
in agreement with their data within 1%. For the 5s5p 
3 P° states, our L-resolved values agree with the CICP 
data at the level of 2%, and our J-resolved result gives 
0 J =1.17xlO 5 for the Sr 5s5p 3 P° state, in consistent with 
the CI+all-order+RPA value 1.05xl0 5 [2l| within 11%. 

C. Black-body Radiation Shift 


Q j ( 3 p f) - Q J ( 3 P 0 °) 

Q J { 3 Pq) [ J 

and 

Q J ( 3po) _ QJfpo) 

Q J ( 3 P$) ’ 1 J 

where Q stands for either a or 7 . Such two fractional 
differences are about 8.9% and 9.6% for a and 28% and 
23% for 7 in the case of In + . Such variations of a J for 
different J components for Sr, as reflected by Eqs. (flOl) 
and ( fTTl) . are about 3.9% using the values in Refs. [ 2 l| 
and [13], and 10.5% using our values. The above values 
of the fractional differences of 5s5p 3 P° and 3 P£ with 
respect to 5s5p 3 Pg are substantial that reveal the im¬ 
portant role of the relativistic effect on a and 7 for the 
states of 5s5p 3 P° of both In + and Sr. 


B. Quadrupole moment and polarizability 

Table IVIII presents 9 and an for the 5s 2 1 5 and 
5s5p 3 P° states of In + and Sr. The L-resolved 
9 l and a% are obtained by the scalar relativis¬ 
tic calculations (corelO)SD(2in4)SDTQ for In + and 
(core8)SD(2in4)SDTQ for Sr, and the J-resolved 9 J 
and «2 are obtained by the relativistic CC calculations 
(corelO)SD for In + and (core 8 )SD for Sr. The data 
for the 5s5p 3 P° state are obtained by the relativis¬ 
tic Cl calculations (corelO)SD(2in4)SD<2 for In + and 
(core8)SD(2in4)SD<2 for Sr. All calculations are carried 


Finally, our results for the dipole polarizabilities, the 
second dipole polarizabilities, and the quadrupole polar¬ 
izabilities are used to estimate the BBR shift in the clock 
transition frequency of In + and Sr. The BBR shift can 
be written in the form Emi 


SE BBR = ~Aa(E 2 E1 ) - ^(E 2 e1 ) 2 - ^Aan(E 2 E2 ), 

( 12 ) 

where (E e1 ) and (E E2 ) are the averaged electric fields 
induced by the electric dipole El and the electric 
quadrupole E2 and they are respectively 


and 


(El 1 ) 


4ti- 3 c4 k B T 4 
15 1 E h > 


(E 


E ’ll 


87t 5 q 4 fc^r 6 
189 1 E h ’ 


(13) 


(14) 


In the above, af s is the fine structure constant, 
ksT/Eh « 10 -9 for the temperature T=300 K, k B is the 
Boltzmann constant, Eh is the Hartree energy, and Aa, 
A «2 and A 7 , expressed in atomic units, are respectively 
the differences of the dipole polarizability, quadrupole 
polarizability and dipole hyperpolarizability between 5s 2 
1 Sq and 5s5p 3 Pq of In + and Sr. In Eq. (TH?1) we have 
neglected the dynamic fractional correction to the total 
shift [13 and assume that the contribution of the hyper¬ 
polarizability to the BBR shift can be approximated by 
the AC-stark shift {E e1 ) 2 for a given electric field. Us¬ 
ing the data of a , 7 , and a-n obtained by the relativistic 
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CC calculations in Tables dun and IyTTI for the states of 
5s 2 1 5 ' 0 and 5s5p 3 P 0 ° of In+ , the BBR shifts due to a, 
« 2 , and 7 are determined to be 0.017, 8.33xlO -10 , and 
1.93xl0 ~ 17 Hz, respectively, for the In + clock transition 
frequency; and using the data of a, 7 , and «2 obtained 
by the relativistic CC calculations in Tables HUE] and 
IVIIl for the states of 5s 2 1 Sq and 5s5p 3 Pg of Sr, the BBR 
shifts due to a, 7 , and ai are determined to be 2.09 and 
5.82 xlO -8 , and 1.69xlO -15 Hz, respectively, for the Sr 
clock transition frequency. It is clear that the contri¬ 
butions from 7 and 0:2 are far less important than that 
from a in the BBR shifts and can thus be safely omitted 
according to the current precision of the quoted 10 -18 
uncertainty of the time frequency standard. 

V. CONCLUSION 

In summary, we have calculated a, 7 , 8 , and a-i for 
the ground state 5s 2 1 5o and the low-lying excited states 
5s5p 3 Pq ) of In + and Sr by using the finite field method. 
A satisfactory accuracy is achieved through convergence 
studies for the basis sets and sufficient inclusion of the 
electron correlations. This method can also be applied 
with similar accuracy to the calculations of the polariz¬ 
abilities of atomic cores, see Supplemental Material at 
[URL will be inserted by publisher] for a, 7 , and a-i of 
In 3+ and Sr 2+ . Thus, it will be useful to employ the fi¬ 
nite field method to perform a fast evaluation of required 
properties, especially when high-precision experimental 
studies and sophistical sum-over-state calculations are 
not available or all available results are not in complete 
agreement. It is noteworthy to mention that the errors of 
the finite field calculations need to be examined carefully 
and minimized for such applications, which requires a de¬ 
tailed knowledge about the rate of convergence of basis 
set and electron correlations for a property of interest. 

In our finite field calculations, we have investigated the 
influences of the basis set and the level of electron corre¬ 
lations on the computed properties. For the case of In + , 
the convergence for a can be easily reached when the ba¬ 
sis set is increased up to A=4£ and the core electrons 
included into the correlation calculation are increased up 
to 3d4s4p4d. The dominant correction to a is from A Pt, 


and thus a more accurate evaluation of the contribution 
of the triple excitation is needed for higher accuracy. For 
example, the relativistic CCSDT calculation with the 4£ 
basis set may give more accurate results, although such 
a study is prohibited at this moment given the present 
computing resources. For the case of Sr, the convergence 
of a is highly sensitive to the quality of the basis set 
because a neutral atom has more diffused electron den¬ 
sity in an external held than a positive charged ion. In 
this case, the predominant factor will be the expansion 
and optimization of the basis sets beyond X=4(. The 
higher-order corrections of the quadruple excitation of 
the electron correlation and the Breit and QED correc¬ 
tions should be roughly two- and one-order of magnitude 
smaller than A Pt in the relativistic CC calculations and 
accordingly the resulted uncertainties should be at the 
level of 0 . 01 % and of 0 . 1 %, which are very small. 

Our fully relativistic calculations have shown that the 
5s5p 3 Pq i 2 states of In + and Sr have obviously differ¬ 
ent values of polarizabilities, which reflects the important 
contributions of the relativistic and spin-orbit coupling 
interactions. Some general trends about the sole effect of 
the spin-orbit coupling are worth noticing through com¬ 
parative studies using the fully and scalar relativistic ap¬ 
proaches. The fractional difference between the L- and 
J-resolved a of 5s5p 3 Pq is about 8% for In + , but only 
2% for Sr, implying that the effect of the spin-orbit in¬ 
teraction on a of In + is stronger than Sr. 
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